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Based on micromagnetic theory we have derived analytical expressions for the magnetic small- 
angle neutron scattering (SANS) cross section of a two-phase particle-matrix-type ferromagnet. The 
approach — valid close to magnetic saturation — provides access to several features of the spin struc- 
ture such as perturbing magnetic anisotropy and magnetostatic fields. Depending on the applied 
magnetic field and on the magnitude Hp of the magnetic anisotropy field relative to the magnitude 
AM of the jump in the longitudinal magnetization at the particle-matrix interface, we observe a 
variety of angular anisotropies in the magnetic SANS cross section. In particular, the model explains 
the "clover-leaf" -shaped angular anisotropy which was previously observed for several nanostruc- 
tured magnetic materials, and it provides access to the magnetic interaction parameters such as the 
average exchange-stiffness constant. It is also shown that the ratio Hp/ IS.M decisively determines 
the asymptotic power-law exponent and the range of spin-misalignment correlations. 

PACS numbers: 61.05.fd, 61.05.fg, 75.25.-j, 75.75.-c 
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I. INTRODUCTION 

Magnetic small-angle neutron scattering (SANS) is one 
of the most important techniques for microstructure de- 
termination in magnetism and magnetic materials re- 
search. The essential advantages of the SANS method 
may be summarized by stating that (i) SANS provides 
nanometer-scale information (^ 1 — 100 nm) from within 
the bulk of a sample and that (ii) SANS experiments can 
be conducted under rather flexible conditions and differ- 
ent sample environments (temperature, electric and mag- 
netic field, pressure, neutron polarization, time-resolved 
data acquisition, etc.). For instance, SANS has previ- 
ously been used for investigating the anisotropy of the 
structure factor and the relaxation dynamics of magnetic 
fluids, '''^■^ the internal structure and morphology of mag- 
netic nanoparticles and nanowire arrays, ^^^ the vortex 
lattice of superconductors, ^° or the spin distribution of 
nanocrystalline^^"^^ and amorphous^^"^^ metals. 

The scattering contrast in magnetic SANS arises 
from deviations of both the magnitude and orientation 
of the local magnetization vector field M(r) from its 
mean value. The phenomenological continuum theory 
of micromagnetics^^"— allows one to compute the mag- 
netic microstructure of a material, i.e., the variation 
of the direction of M as a function of the position r 
and time. Since the elastic magnetic differential scat- 
tering cross section depends on the Fourier coefficients 
of the magnetization, it is, in principle, straightforward 
to use micromagnetic theory for predicting and model- 
ing of magnetic SANSi^"— However, due to the nonlin- 
earity of the governing equations — Brown's equations of 
micromagnetics — it is rather difficult to obtain analytical 
closed- form expressions for the cross section and, indeed, 
such solutions are limited to the approach-to-saturation 



regime where the micromagnetic equations can be lin- 
earized. 

Based on the pioneering work of Kronmiiller et al. 
(Ref. 1321 ) we have previously derived analytical expres- 
sions for the differential magnetic scattering cross section 
of ferromagnets with uniform values of the saturation 
magnetization and the exchange interaction, but with 
a highly nonuniform magnetocrystalline and/or mag- 
netoelastic anisotropy. ■^■^'■^■* Examples for such materials 
are nanocrystalline or cold-worked metals which exhibit 
nanoscale variations in the direction and/or magnitude 
of the magnetic anisotropy, leading to a highly inho- 
mogeneous spin structure along with a strongly field- 
dependent magnetic SANS signal. As summarized in 
Ref. 1351, analysis of experimental SANS data on nanocrys- 
talline cobalt and nickel and on cold-worked coarse- 
grained nickel yields information on the strength and 
spatial structure of the magnetic anisotropy field, on 
the magnitude and field dependence of the magnetostatic 
field, on the exchange constant, and on the characteristic 
length of the spin misalignment. 

It is the purpose of this article to extend the micromag- 
netic scattering approach to materials having nonuni- 
form saturation magnetization. Prominent examples 
are two-phase hard and soft magnetic nanocomposites, 
which are used e.g. in electronics devices, transformers 
or motorsi^^*2I Here, the jump in the magnitude of the 
magnetization at the particle-matrix interface gives rise 
to a magnetostatic stray field which represents a domi- 
nating source of spin disorder ;3S, 

The paper is organized as follows: To start with, we 
summarize in Sec. |TT] well-known results for the mag- 
netic SANS cross section of a fully saturated two-phase 
particle-matrix- type ferromagnet. Section IIIII describes 
the micromagnetic model and provides approximate ex- 
pressions for the Fourier coefficients of the magnetization, 
which are used in Sec. |lV]in order to compute the mag- 
netic SANS cross section; here, we focus on the two most 



relevant scattering geometries where the externally ap- 
plied magnetic field is perpendicular or parallel to the 
incident neutron beam. Section |V] discusses the results 
in real space in terms of the corresponding spin-spin cor- 
relation function. Finally, Sec. IVII summarizes the main 
findings of this work. 



II. SATURATED STATE 

It is our aim to derive an approximate expression for 
the magnetic microstructure and the ensuing magnetic 
SANS cross section of a two-phase nanocomposite ferro- 
magnet, i.e., a bulk material which consists of a distri- 
bution of ferromagnetic nanoparticles that are embed- 
ded in and magnetically exchange-coupled to a ferro- 
magnetic matrix; all particles are assumed to possess 
a saturation magnetization of Mp^ whereas the satura- 
tion magnetization of the matrix phase is denoted with 
Mm- We remind the reader that for small-angle scat- 
tering the discrete atomic structure of both particle and 
matrix phase is generally of no relevance. Therefore, as 
far as magnetic SANS is concerned, the magnetization 
state may be represented by a continuous magnetization 
vector field M(r). Before discussing in detail in Sec. IIIII 
the micromagnetic model of such a two-phase nanocom- 
posite structure, we find it instructive to consider first 
the (somewhat idealized) fully saturated state. 

If, for instance, the external magnetic field Hg || e^ 
is large enough to fully align the magnetic moments of 
the sample, then the local magnetization simplifies to 
M(r) = (0,0, Mz), where the longitudinal component 
Mz is cxplicitely dependent on the position r inside the 
material, i.e., 

M,(r) = Mp u{v) + M„ [1 - u{v)] , (1) 

where u(r) = 1 inside the particle phase and u(r) = 
inside the matrix phase. By assuming that the function 
Mz{r) can be expressed as a Fourier integral, 

Mz{r) = —^ I MzicO exp(zqr) d\ , (2) 



(27r)3/2 



M.(q) = 



1 



Mz{r) exp(-iqr)(i^ 



(3) 



(27r)3/2 

it is straightforward to obtain the following limiting 
cases: At q = 0, one obtains the macroscopic saturation 
magnetization Mg of the sample (which can be measured 
with a magnetometer). 



Ms^{Mz) ^V-^ / Mzir)d^r 



(2^) 



3/2 , 



^ -Mz (q = 0) = Mp p + Mm (l-p) , (4) 

where V is the sample volume and p equals the volume 
fraction of the particles. At q 7^ 0, we obtain^S 



*^-(q)= (^^:j37i^p(q) 



(5) 



where AM = Mp — Mm denotes the jump in the magne- 
tization magnitude at the particle-matrix interface and 
the function /p(q) contains information about the size, 
the size distribution, the shape and the arrangement of 
the particles (compare Eq. p2|) below). 

In the monodisperse and dilute limit we find for a single 
spherical particle (with radius R) 



Mz{q) 



^M nm 



(27r)3/2 



qR 



(6) 



where Vp = ^R^ and ji{qR) denotes the spherical Bessel 
function of first order. By using Eq. ([B]) we can immedi- 
ately write down (in the monodisperse and dilute limit) 
the magnetic SANS cross section dY.M /dVl for a collec- 
tion of Np saturated particles in a saturated matrix. For 
the particular scattering geometry where the wave vec- 
tor ko of the incident neutron beam is perpendicular to 
the applied magnetic field Hq, we obtain the well-known 
expression 



^(q) = ^ApLg^/f^(9i?)sin^^, 



(7) 



where Ap^^ = 6|^(AM)2 represents the magnetic 



scattering-length density contrast, F{qR) = 3 """„'' is 



ji(qR) 



the form factor of the sphere, and 9 denotes the angle be- 
tween the scattering vector q and Hg (compare Eq. P^ 
below) . 

Deviations from the saturated state result in the emer- 
gence of transversal magnetization components. The as- 
sociated so-called spin-misalignment scattering will be 
addressed in the following sections. 



III. MICROMAGNETIC MODEL 

Analytical micromagnetic calculations of the type pre- 
sented here have already been carried out by other au- 
thors (e.g., Refs. I39l - l4l[) . We employ the micromagnetic 
approach for computing (in the high-field limit) the dif- 
ferential magnetic scattering cross section (see Sec. lIVp . 
In the following we summarize the basic micromagnetic 
equations. 

We are interested in the elastic spin-misalignment scat- 
tering which results from the static magnetic microstruc- 
ture of a two-phase nanocomposite sample. Therefore, 
we start our analysis by writing down Brown's balance- 
of-torques equation^^"— 



M(r) X Heff(r) = 



(8) 



which expresses the fact that at static equilibrium the 
torque on the magnetization vector field ]V[(r) due to an 
effective magnetic field Hoff(r) vanishes everywhere. In 
the micromagnetic model we assume a uniform exchange 
interaction but an explicitely wave- vector-dependent lon- 
gitudinal magnetization (see below). The effective field 



Heff(r)=Ho + Hd(r)+Hp(r) + 



2A 



V2M(r) (9) 



is composed of a uniform applied magnetic field Hq, of 
the magnetostatic field Hc;(r), of the magnetic anisotropy 
field Hp(r), and of the exchange field [last term on the 
right hand side of Eq. (0)]; A^o = 47r x 10"'^Tm/A is 
the permeability of free space, V = e^^ d/dx + ey d/dy + 
Gz d/dz^ where e^;, ey, and Gz represent the unit vectors 
along the Cartesian laboratory axes. The parameter A 
denotes the exchange-stiffness constant. 

In the following we assume the material to be nearly 
saturated along Hq |1 e^, i.e., we write 



M(r) = M^(r) e^, + My{v) By + Mz{v) e^ 



(10) 



with Mx ^ Mz and My ^ Mz (small-misalignment ap- 
proximation). The local saturation magnetization is as- 
sumed to differ only slightly from its spatial average, i.e., 
Mz{y) = {Mz) = Ms. Note that the jump l\M in the 
longitudinal magnetization enters the calculation in q- 
space via the expression for the Fourier coefficient of the 
magnetostatic field (see below) . Furthermore, we assume 
that the anisotropy-energy density uj = a;(r, M) depends 
only linearly on the components of the magnetization)^ 
As a consequence, the resulting anisotropy field Hp = 
~fi^^{duj/dMx, du}/dMy, dw/dMz) is independent of M 
and, therefore, also independent of the applied magnetic 
field, implying that near saturation Hp — Hp(r). Due to 
the micromagnetic constraint |M| = Mg, an anisotropy- 
energy density of the form w = u]{r,Mx,My,Mz) may 
be re-expressed as cj = uj{r,Mx,My) with the conse- 
quence that only two independent components of Up ex- 
ist. In the approach-to-saturation regime, when M is 
nearly aligned parallel to the external magnetic field Hq , 
only those components of Up which are normal to Hq are 
physically effective in producing a torque on the magne- 
tization. 

Basic magnetostatics prescribes that V • (Hq + Hd) = 
—V • M and that V x (Ho + Hd) = 0. The magnetostatic 



field H(;(r) can be written as the sum of the surface de- 
magnetizing field H^ and of the magnetostatic field H^ 
which is related to volume charges, i.e., H^ = H^ + H^. 
In the high-field limit and for samples with an ellipsoidal 
shape with Hg directed along a principal axis of the ellip- 
soid, one may approximate the demagnetizing field due to 
the surface charges by the uniform field H^ = —NMgez, 
where N denotes the corresponding demagnetizing fac- 
tor. In Fourier space (at q 7^ 0) the above magnetostatic 
relations suggest the following expression for the Fourier 
coefficient h^(q) of H;j(r)^ 



h:^(q) 



q[qM(q)] 



(11) 



H^(r) 



(27r)3/2 



h^(q)exp(zqr)d3g. (12) 



M(q) represents the Fourier transform of the magnetiza- 
tion M(r), 

M(r) = ^^-^ J M(q) exp(iqr) d^q . (13) 

Likewise, the Fourier transform h(q) = 

(/ia;(q), ft,j,(q), 0) of the magnetic anisotropy field 
Hp(r) is introduced as 



Hp(r) = 3^^ / h(q) exp{iqr)d^q 



(14) 



The details of the sample's microstructure (e.g., grain 
size, lattice strain, crystallographic texture) are included 



in Hp(r).34 



By inserting Eqs. ([TO]) - ([Ml) into Eqs. (g]) and ^, 
we obtain for a general orientation of the wave vector 
q = {qx,qy,qz) the following expressions for Mx{q) and 
M,(q): 



Mxiq) = M, 



My{q) = M, 



{hx-Mz^-^){ 



Hes + Ms-p 



)-Af,^( 



Af,^4^ h„-M,^ 



{hy - Mz^) (ffcff + Afs J) -Ms^ (/i. 



Mz^ 

z q2 



Hrl 



yHcs 






(15) 
(16) 



r 



The terms in Eqs. P^ and ((TE| which contain the 
Fourier coefficient Mz{q) oc AM model the influence of 
the two-phase magnetic microstructure on the magnetic 
SANS and are not contained in the corresponding ex- 
pressions for the single-phase case (compare Eqs. (2.15) 
in Ref . |35I ) . The quantity 



2 „2\ 



H,s{q,H,)^H,{l + lj,q 



(17) 



is the effective magnetic field, which depends on the in- 
ternal field Hi = Hq — NMs, on q = |q|, and on the 
exchange length of the field IniHi) = y^2A/{fioMsHi); 
the length scale Ih characterizes the range over which 
perturbations in M decay (see Sec. |V] below) .—i^ 



In the derivation of Eqs. ([T5|) and (|T6l) terms of higher 
than linear order in Alx or My have been neglected, in- 



eluding terms such as H^ ^ M^ or H^ ■ My where i G 
{a;,y, z}.— Besides this small-misahgnment approxima- 
tion, we have introduced two further approximations: (i) 
The exchange interaction is assumed to be homogeneous, 
i.e., jumps in A between the two magnetic phases have 
not been taken into account. Such an approximation 
is permissable as long as exchange fluctuations are not 
too large, in particular, for soft magnetic materials.-^ 
(ii) The function ^/^(q), which models the influence of 
the two-phase magnetic microstructure on the magnetic 
SANS, is introduced into our theory only in q-space 
via h^(q) = -q[qM(q)]/g2 [Eq. (dJ)]. This is an ap- 
proximation, since in real space we have assumed that 
Mz = Ms = constant, and hence Mziq) ex 6{q) would 
result, as is appropriate for a homogeneous single-phase 
material. However, by explicitly considering the q 7^ 
Fourier coefficients of Mz [Eq. ([5])], it becomes possible 
to straightforwardly include the jump in the longitudinal 
magnetization at the particle-matrix interface, and one 
avoids the otherwise necessary calculation of convolution 
products 1^"— 

For the following discussion of magnetic SANS it is 
of interest to consider special projections of Eqs. P^ 
and ([TO)) into the plane of the 2D detector. The two 
scattering geometries which are of particular relevance 
to experiment have the external magnetic field Hq either 
perpendicular or parallel to the wave vector ko of the 
incoming neutron beam. For ko || e^^ and Hq || e^, the 
scattering vector can be approximated as q = (0, g^, q^), 
i.e., Qx — 0, and Eqs. (IT5t and P^ reduce to 



M,(q) = M, 



hxjq} 



~ /iJq)-M^(q)2^ 
My{ci)^Ms^^ ^^^^^^ 



H,s + Ms^ 



(18) 
(19) 



Equations ((20)) and (|2T|) are similar to Eqs. (2.15) in 
Ref. m. 

The high- field solutions for Mxiq) and My{q) can 
be seen as a sum of products of components of the 
anisotropy-field Fourier coefficient h(q) and Mz{q) with 
micromagnetic functions which contain the effective mag- 
netic field HcS and terms that depend on the orientation 
of the wavevector (e.g., Ms ql/q^)- The convolution the- 
orem then implies that the magnetic microstructure in 
real space, M(r), is tantamount to a complicated con- 
volution product between the corresponding real-space 
functions. As a consequence, sharp features in the nu- 
clear or anisotropy-field microstructure are washed out 
and smoothly-varying magnetization profiles are at the 
origin of the related spin- misalignment scattering (com- 
pare, e.g.. Figs. 2 — 4 in Ref. Hj). Consistent with this 
notion is the observation of power-law exponents signif- 
icantly larger than 4 (see Fig. 2] below) and the finding 
that the slope of the correlation function at the origin 
vanishes (see Fig. [S] below) 4^ 



IV. SANS CROSS SECTION FOR 
UNPOLARIZED NEUTRONS 



Although the following derivation of the magnetic 
SANS cross section is aimed to be self-contained, we refer 
the reader to Refs. w3 and wj for further details. 



For ko II Ho || e^, q = {qx,qy,0), i.e., qz = 0, and the 
results for the Fourier coefficients simplify to 



_ hx (i?eff + Ms$)-hy Ms ^ 

M.(q) = Ms -^ ^ ,.+,.x , (20) 

H,s [H,s + Ms'-^ 

_ hy (h,s + Ms^)- hx Ms^ 

My{<l) = Ms -^ J ' ^ ^,,„,, . (21) 



H, 



off 



(^cff 



M< 



Qi+yi 



J 



ko -L Ho 



For the transversal scattering geometry, one can ex- 
press the total nuclear and magnetic SANS cross section 
dJ^/dil at scattering vector q as^ 







+ |Mj/|2cos^6i + |Af, 



Psin^i 



{MyM; + M*Mz) sin 9 cos 



(22) 



r 



Here, V is the scattering volume, 6^ = 2.9x 10® A ^m ^, 
A^(q) is the nuclear scattering amplitude, and 9 denotes 
the angle between Ho and q = g (O,sin0, cos6'); c* is a 
quantity complex-conjugated to c. The magnetic form 



factor in the expression for the magnetic scattering am- 
plitude (ex hu) was set to unity, which is permissible 
along the forward direction. 



By inserting Eqs. p8|) and ([T9| into Eq. ([22]) one can 
express (E/dVt as 






where 



A^l 






lAfJ^sin^ 



(23) 



(24) 



if 



represents the (nuclear and magnetic) residual SANS 
cross section, which may be measured at complete mag- 
netic saturation. Subtraction of dYi-ccs/d^ from a mea- 
surement of dTi/dCl at a lower field yields the spin- 
misalignment SANS cross section dT,M/dfl, which con- 
tains the scattering contributions due to the misaligned 
spins, 



dE 



M 



dfl 



{q) = Sniq) RHiq,0,H,) + SM{q) RMiq,9,Hi 



_ (25) 

Following Weissmiiller et al. (Ref. |3J), we have here 
introduced the scattering function of the anisotropy field. 



SH{<i)^—bj,h'{<l), 



(26) 



the scattering function of the longitudinal magnetization, 

87r3 



5A./(q) = ^6?,M2(q) 



(27) 



and the corresponding (dimensionless) micromagnetic re- 
sponse functions. 



and contain, respectively, information on the strength 
and on the spatial structure of the magnetic anisotropy 
field^ as well as on the magnitude of the jump of the 
magnetization at the particle-matrix interface. The ge- 
ometry of the microstructure is contained in Sh and Sm- 
Both response functions depend on the magnitude and 
orientation of the scattering vector, on the applied field, 
and on the magnetic interaction parameters. 

As mentioned previously, the effects of crystallographic 
texture or of other forms of anisotropy on dT,M/dfl en- 
ter the theory mainly through the magnetic anisotropy 
field Hp(r). In deriving Eqs. (US])"® we have as- 
sumed that the corresponding anisotropy-field Fourier 
coefficient h(q) = (/i(q) cos/3, /i(q) sin/3, 0) is isotropi- 
cally distributed in the plane perpendicular to Hq, in 
other words, the vector h(q) takes on all orientations 
(i.e., angles /3) with equal probability. This assumption 
allows us to average the response functions over /3, i.e., 
l/(27r)/Q ''(...) d/3, which results in Eqs. ^ and (|29l) : 

note that interference terms oc hy M^ = h sin /3 M^ van- 
ish in this averaging procedure (see also Sec. lIVCJ belowV 
The case of a texture in the orientation of the anisotropy 
field is treated in the Appendix of Ref. |3J. 

Furthermore, by assuming that both Fourier coeffi- 
cients h? and M| depend only on the magnitude q of 
the scattering vector q, and not on its orientation 9, one 
may average Eq. ((25|) [and Eq. (l30l) ] with respect to the 
angle 0, i.e., l/(27r) /^ (...) d9, which results in 



dE 



M 



dVL 
where 



(q) - SH{q) RH{q.H,) + Sniq) RM{q,H,) (31) 



RH{q.e,Hi) 



P 



cos^9 



[I +psin 



(28) 



RH{q,Hi)^^(2 



4 V V^+^J ' 



(32) 



„ , „ ^, , B^ sin^ 6' cos** 2p sin^ 6' cos^ 6* ,„„, 

RM{q, e, H{) = 5 — ^ + \ . , (29) 

{l+psm^e) 1+psm 6 



RM{q,Hi 



VT+P-l 



(33) 



where p{q,Hi) = Ms/Hcff{q,Hi). We note that the term 
Sm X Rm in Eq. ([25| [and Eq. (|3T|) ] is not contained in 
the expression for the single-phase material case (homo- 
geneous saturation magnetization), where^ 



Equation (I3ip may be compared to experimental data for 
dYiM /d^ in order to obtain the functions Sh and Sm and 
a value for the average exchange-stiffness constant A. 



dY. 



M 



dn 



{q)^SHiq)RHiq,e,H,) 



(30) 



Near magnetic saturation, Sh and Sm are both ap- 
proximately independent of the applied magnetic field 

I 



B. ko II Ho 



For the longitudinal SANS geometry, dE/dil reads 



dE 
dh 



(q) 



V 



■hi 






\MA^sa?B 



[Mypcos^ 



B + |/V4|2 - [M^M* -f MlMy) sin6'cos6' 



(34) 



where q = (7 (cos 61, sin 0,0) and 9 — Z(e2.,q). Inserting 
Eqs. ((20)) and (|2T|) into Eq. (p4)) and averaging over the 
orientations of the magnetic anisotropy field (angle (3) 
results in 






dYj-r 



-(q) 



dS 



M 



(q), 



d^ ^^' dn 

where the residual SANS cross section equals 






^ 
^^1. 



\M, 



(35) 



(36) 



and the spin-misalignment SANS cross section is ex- 
pressed as 



dS 



M 



dn 



{q) = SHiq)RH{q,H,) , 



(37) 



with an isotropic (i.e., ^-independent) response function. 



RH{q,H^)^ 



(38) 



We note that in the longitudinal SANS geometry 
d'^M /dfl is independent of Mz and equals the expression 
for the single-phase case [compare Eq. (33) in Ref. l33| . 
In other words, the two-phase nature of the underlying 
microstructure is (for kg |1 Hq) only contained in the 
residual SANS, and not in dT,M/dVL. 



C. Comment on the average of dY^M /dQ, over many 
statistically uncorrelated defects 

The sample volume which is probed by the neutrons 
typically contains many defects (e.g., particles), each 
one having a different orientation and/or magnitude of 
the magnetic anisotropy field. For ferromagnets with 
a nonuniform saturation magnetization, each particle or 
crystallite "j" is additionally characterized by the form- 
factor function Mzj. 

In order to discuss the statistical average over the spin- 
misalignment SANS cross sections of many defects, we 
ignore for the moment the AM fluctuations, i.e., we 
consider the case of a homogeneous single-phase ferro- 
magnet with uniform values of A and Ms- Here, the 
dominating source of spin disorder (in the approach-to- 
saturation regime) is related to spatially inhomogeneous 
magnetic anisotropy fields. We follow the arguments of 
Weissmiiller— 1^ and assume that the total anisotropy- 
field Fourier coefficient of the sample, h(q), can be ex- 
pressed as the sum of the anisotropy-field amplitudes of 
the individual defects. 



h(q) = ^h,(q) 



(39) 



If the hj of the individual defects are statistically uncor- 
related (random anisotropy), then the expectation value 



of |h(q)p reduces to the sum over the individual |hjp, 
|h(q)p=^|h,(q)r (40) 



Since — within the linear micromagnetic approximation — 
dYj^i/dO, is proportional to |h(q)p [compare Eq. [30], it 
immediately follows that the above additivity also trans- 
fers to the total spin-misalignment SANS cross section of 
the sample. In other words. 



dS 



M 



dn 



E 



dS 



M-3 



dfl 



22 ^H,j Rh,j 



(41) 



when homogeneous single-phase ferromagnets are 
considered. •^■^'•^■^ This is in contrast to nuclear SANS 
or SAXS, where the decomposition of the overall cross 
section into the sum over cross sections of individual 
particles is only permissible for small particle volume 
fractions. 

Including now the AM fluctuations into the discus- 
sion, inspection of Eq. (|T9t shows that interference terms 
oc hyj Mzj appear in the spin-misalignment SANS cross 



section [compare Eq. (1221) ]. However, in deriving the fi- 
nal expression for dYjM/dVt [Eq. (P5|) ] these terms cancel 
when the averaging procedure over the (isotropic) orien- 
tations of the anisotropy field (angle (i) are carried out, 
with the result that the scattering contributions due to 
inhomogeneous anisotropy fields Sh x Rh and magneto- 
static fluctuations Sm x Rm are additive. Taking into 
account intcrparticle interference effects, the general ex- 
pression for \Mz\'^ (at q 7^ 0) may be expressed as^ 



|M,(q) 



{AMf 
87r3 



E^PJ^j(^) exp(-iqrj 



(42) 

where Np, Vpj, Fj and r-,- denote, respectively, the num- 
ber of particles, the volume, the form factor and the posi- 
tion vector of particle j [compare Eq. ^]. Equation (|i^ 
can be employed in Eqs. ([25)) or ([31)) in order to analyze 
experimental data. 



D. Graphical representation of dT,M /dO. 



In order to graphically display dY^M/dO. [Eq. ([^ ] it is 
necessary to specify particular models for the Fourier co- 
efficients h{q) and Mz{q). For Mz{q), we use for simplic- 
ity the form factor of the sphere [Eq. ([5])] , which implies 
the neglect of interparticle interactions (dilute limit). We 
note, however, that the inclusion of the structure factor of 
the material, of other particle shapes (form factors), or of 
the particle-size distribution is straightforward [compare 
Eq. ([1^ ].^- For h{q), we employ the model introduced by 
Weissmiiller et al. (Eq. (43) in Ref. |33I) . which considers 
a nanocrystalline material composed of spherical parti- 
cles with a constant magnitude but random orientation 
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FIG. 1. Micromagnetic response functions Rh (Eq. (|32[l : 
solid line) and Rm (Eq. (|33|) : dashed line) versus p{q,Hi) = 
Ms/HcB{q,Hi) (log-log scale). 



of Hp. The resulting expression for h(q) is the form fac- 
tor of the sphere, Eq. (jH]), where AM is replaced by the 
magnitude Hp of the magnetic anisotropy field. For the 
particle radius, we assume a value of i? — 5 nm for both 
Mz{qR) and h{qR). Since, then, Mz{qR) and h(qR) and, 
hence, Sm and Sh are equal except for the pref actors, 
the ratio Hp/AM determines the properties of dYjM /<Kl. 
Unless otherwise stated, we have used the following ma- 
terials parameters: A = 2.5 x 10"-^^ J/m; ^iqMs = 1.5 T; 
/ioAM = 0.25T. 

The spin-misalignment SANS cross section for ko -L 
Hq [Eq. (|25l) ] contains scattering contributions due to the 
magnetic anisotropy field, Sh x Rh, and due to jumps 
in the longitudinal magnetization, Sm x Rm- In Fig. [1] 
we plot both response functions Rr [Eq. ([5^ ] and Rm 
[Eq. (|33p ] as a function of the dimensionless parameter p. 
Assuming that Sr and Sm are of comparable magnitude, 
it is seen that at large applied fields or large momentum 
transfers (when p <C 1), dT^M/dfl is dominated by the 
term Sm x Rmi whereas at small fields and small mo- 
mentum transfers (when p 3> 1), d'EM/d^ is governed 
by the Sh x Rh contribution (see Fig. [5] below). 

Figure O qualitatively displays the applied-field de- 
pendence of dTiM/d^ for ko J- Hq [Eq. ((25|) ] and for 
a ratio of Hp/AM = 1 [Sh — Sm)] the results are 
compared with the dY^M/dU, of a single-phase mate- 
rial with a uniform saturation magnetization [Eq. pop ]. 
The dYjM/dO. for the two-phase case [Figs.^a) — (d)] re- 
veal a strongly field-dependent angular anisotropy. At 
the largest fields, the pattern exhibits maxima roughly 
along the diagonals of the detector — the so-called "clover- 
leaf" anisotropy — previously observed in the Fe-based 
two-phase alloy NANOPERM (compare, e.g., Fig. 3 in 
Rcf. 38). The position of the maxima in dY^M/dU, de- 
pend on q and Hi (see also Fig. 11 in Ref. [3l|). Such a 
type of clover-leaf anisotropy cannot be reproduced by 
the dT,M/dn for the single-phase case [Figs. [IJe) — (h)]. 
Here, at large q and Hi, one observes an elongation of 
the spin-misalignment scattering along the field direction 



[due to the cos^ 9 term in Eq. ([28|) ]. with a "flying-saucer- 
type" pattern taking over at small q and Hi [due to the 
sin^ 6 term in the denominator of Eq. ([^5]) ]. The sharp 
spike in Figs.UJa) and (e) is due to the magnetostatic in- 
teraction and was first predicted by Weissmiiller et aL^^ 

Figure [3] shows dT,M/dn for ko -L Hq [Eq. 1^] at 
a fixed internal field of 2.0 T but for different magni- 
tudes of the magnetic anisotropy field Hp relative to 
the jump AM in the magnetization magnitude at the 
particle-matrix interface. One can clearly observe a 
transition from dipole-field dominated scattering, with 
a characteristic clover- leaf- type pattern [Fig. [3I^a)], to 
a more anisotropy-field dominated cos^ 9-type angular 
anisotropy of dT,M/dQ [Fig. ^c)]. 

Regarding the asymptotic power-law dependence of 
dYjM /d^: For particles with sharp interfaces, both h?{q) 
and M^{q) vary asymptotically as q^"^, as does the func- 
tion -ffpff^ [compare Eq. p?|) ]. Taking these dependencies 
into account, it is readily verified that the anisotropy- 
field contribution to dY^M / dO, varies as Sh x Rh oa g^^, 
whereas Sm x Rm oc q^^. Therefore, depending on the 
relative magnitude of both contributions to d^M /dfl, one 
observes different asymptotic power-law exponents (see 
Fig. 2]) . We note that other models for the anisotropy- 
field microstructure may result in different exponents. 



V. CORRELATION FUNCTION OF THE SPIN 
MISALIGNMENT 

The results for dYM /dVt [Eq. (|3ip ] can be used in order 
to compute the autocorrelation function C{r) of the spin 
misalignment, according tc4^i^i^ 



C{r) ex - 



dY 



M 



dfl 



(q) sm{qr) q dq . 



(43) 



Equation (|43|) has been solved numerically and the re- 
sults for C{r) at ^oHi = LOT and for different ratios 
of Hp/AM are shown in Fig. [5l The corresponding field 
dependence of the so-called correlation length Ic of the 
spin misalignment, which quantifies the size of inhomo- 
geneously magnetized regions, can be seen in Fig. [S] 

In agreement with the previous studies (e.g., Refs. l43l 
and |4J), we find strongly field-dependent spin- 
misalignment correlations. The slope of the correla- 
tion function at the origin vanishes, which is in line 
with the absence of a sharp interface in the magnetic 
microstructure.— 1^ The ratio Hp/AM decisively deter- 
mines the characteristic decay length Ic'- Increasing 
Hp/AM results in the emergence of more long-range 
magnetization inhomogeneities, whereas AM dominated 
perturbations in the spin structure decay on a smaller 
length scale. For Hp/AM ^ 1, the relation lc{Hi) can 
be described by the phenomenological expression (solid 
line in Fig. [HI'*" 
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FIG. 2. (Color online) (upper row) Contour plots of dTiM/dD. (in arbitrary units) for ko -L Ho [Eq. (I25p ]. Ho is horizontal. 
Hp/AM ^ 1 {Sh = Sm). Hi values (in T) from (a) to (d): 0.01; 0.2; 1.0; 10.0. (lower row) Field dependence of dY.M/dQ. for 
ko 1- Ho for a homogeneous single-phase ferromagnet [Eq. (|30p ]. where AAf = 0. Field values and materials parameters in (e) 
to (h) are the same as in (a) to (d). Red color corresponds to "high intensity" and blue color to "low intensity". 
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FIG. 3. (Color online) Crossover from magnetostatic to 
anisotropy-field dominated scattering. Contour plots of 
dTiM/dO, (in arbitrary units) for ko _L Ho [Eq. (|25|) ] at a 
fixed internal magnetic field of ^oHi = 2.0 T. Ho is hori- 
zontal. Values of Hp/AM from (a) to (c): 0.2; 1.6; 8. Red 
color corresponds to "high intensity" and blue color to "low 
intensity" . 



Equation (|44|) embodies the convolution relationship 
between the anisotropy-field microstructure (through 
the field- independent parameter R) and micromagnetic 
functions (through the field-dependent exchange length 
y^2A/{iioMsHi). Irrespective of the value of Hp/AM, 
it is seen that at large fields Ic approaches the particle 
radius (dashed line in Fig. |6|). 



VI. SUMMARY AND CONCLUSIONS 

Using the continuum theory of micromagnetics we have 
derived in the approach-to-saturation regime analytical 
expressions for the magnetic small-angle neutron scatter- 
ing cross section of a two-phase particle-matrix-type fer- 
romagnet. For the particular scattering geometry where 
the applied magnetic field is perpendicular to the incom- 
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FIG. 4. Variation of the power-law exponent n in dY^M/dO. = 
K/q" with the ratio Hp/AM {^oHi = LOT) (ko ± Ho). The 
fits of the above function to the simulated data [Eq. (|3ip ] have 
been restricted to the interval 1.0 nm~ < q < 2.0 nm~ . Line 



is guide to the eye. 



ing neutron beam, the results for the spin-misalignment 
cross section dY^M /d^ [Eq- ([25]) ] exhibit a variety of angu- 
lar anisotropics that are fundamentally different from the 
conventional sin or cos^ 0-type patterns. In particular, 
by explicitly taking into account the wave-vector depen- 
dence of the longitudinal magnetization, novel terms ap- 
pear in dTiM/dn which give rise to maxima roughly along 
the diagonals of the detector ( "clover-leaf" anisotropy) , 
in agreement with experiment. Besides the value of 
the applied magnetic field, it is the ratio of the mag- 
netic anisotropy field Hp to the jump AM in the Ion- 
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FIG. 5. (Color online) Normalized correlation function C(r) 
of the spin misalignment at yLoHi — LOT and for different 
ratios of Hp/AM (decreasing from top to bottom, see inset) 
(ko -L Ho). Dashed horizontal line: C{r — Ic) ~ exp(— 1). 



gitudinal magnetization at internal interfaces which de- 
termines the properties of dTtM/dfl, for instance, the 
asymptotic power-law exponent, the angular anisotropy, 
or the characteristic decay length of spin-misalignment 
fluctuations. Therefore, analysis of dSjvf/dfi will pro- 
vide information on the strength and spatial structure 
of the magnetic anisotropy field, on AM, and on the 
effective exchange-stiffness constant. Computation of 
the spin-spin correlation function corroborates the exis- 
tence of long-range spin-misalignment fluctuations, with 
a characteristic field-dependent length that in the high- 
anisotropy-ficld limit {Hp/AM 3> 1) can be described by 
Eq. ([44)1 . By contrast, perturbations in the spin struc- 
ture that are due to magnetostatic fluctuations decay on 
a relatively short length scale. 
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